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Measurement of the Decay Constant f D + using Df — > i + i 
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We measure the decay constant f D + using the — > l + v channel, where the £ + designates either 
a n + or a t + , when the r + — » n + V. Using both measurements we find f D + = 274 ± 13 ± 7 MeV. 
Combining with our previous determination of fo+, we compute the ratio f D +/fo+ ~ 1.23 ± 0.11 ± 
0.04. We compare with theoretical estimates. 

PACS numbers: 13.20.Fc, 13.66.Bc 



To extract precise information on the size of CKM 
matrix elements from Bd and B s mixing measure- 
ments the ratio of "decay constants," that are re- 
lated to the heavy and light quark wave-function 



overlap at zero separation, must be well known [l[. 
Recent measurement of 5° mixing by CDF has 
shown the urgent need for precise numbers. De- 
cay constants have been calculated for both B and 
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D mesons using several methods, including lattice 
QCD Here we present the most precise measure- 
ment to date of f D + , and combined with our previ- 
ous determination of f D + @, @|, we find f D +/fo+- 

In the Standard Model (SM) purely leptonic D s 
decay proceeds via annihilation through a virtual 
W + . The decay rate is given by [6[ 
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where M D + is the D+ mass, mg is the lepton mass, 
Gf is the Fermi constant, and |T4J is a CKM matrix 
element with a value of 0.9738 pf. 

In this Letter we report measurements of both 



~v) and B(D+ — > t + v\ when t h 



B(D+ - M _ 

(Z?J — > tt + uv). More details are given in a compan- 
ion paper % The ratio T{D+ -> T+z/)/r(D+ -> 
predicted in the SM via Eq. [1] depends only 
on well-known masses, and equals 9.72; any devia- 
tion would be a manifestation of new physics as it 
would violate lepton universality New physics 
can also affect the expected widths; any undiscov- 
ered charged bosons would interfere with the SM 

W+ 0. 

The CLEO-c detector ll| is equipped to mea- 
sure the momenta of charged particles, identify 
them using dE/dx and Cherenkov imaging (RICH) 
(l2l |. detect photons and determine their directions 
and energies. We use 314 pb _1 of data produced 
in e + e~ collisions using CESR near 4.170 GeV. 
Here the cross-section for our analyzed sample, 
D* + Dj +Df D*~ , is ~1 nb. Other charm produc- 
tion totals ^7 nb 13], and the underlying light- 
quark "continuum" is ~12 nb. We fully reconstruct 
one D~ as a "tag," and examine the properties of the 
Df. (Charge conjugate decays are used.) Track se- 
lection, particle identification, 7r°, ry, and K% criteria 
are the same as those described in Ref. [4j, except 
that RICH identification now requires a minimum 
momentum of 700 MeV/c. 

Tag modes are listed in Table HI For resonance de- 
cays we select intervals in invariant mass within ±10 
MeV of the known mass for rf — » 7r + 7r~?7, ±10 MeV 
for <j> -> K+K~, ±100 MeV for K* a -> R--K+, and 
±150 MeV for p~ — > ir~ir°. We require tags to have 
momentum consistent with coming from D S D* pro- 
duction. The distribution for the K + K~ir~ mode 
(44% of all the tags) is shown in Fig. [TJ 

To select tags, we first fit the invariant mass dis- 
tributions to the sum of two Gaussians centered 
at Md s . The r.m.s. resolution (<j) is defined as 
a = /i<7i + (1 — /i)<72, where <T\ and a 2 are the in- 
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FIG. 1: Invariant mass of K + K~tt~ candidates after 
requiring the total energy to be consistent with the beam 
energy. The curve shows a fit to a two-Gaussian signal 
function plus a polynomial background. 



TABLE I: Tagging modes and numbers of signal and 
background events, within cuts, from two-Gaussian fits 
to the invariant mass plots, and the number of 7 tags in 
each mode, within ±2.5<r from a fit to the signal Crys- 
tal Ball function (see text) and a 5th order Chebychev 
background polynomial and the associated background. 



Mode 


Invariant Mass 


MM* 






Signal 


Bkgrnd 


Signal 


Bkgrnd 


K+K'ty' 


13871±262 


10850 


8053± 211 


13538 


K%K~ 


3122L79 


1609 


1933L88 


2224 




1609 ± 112 


4666 


1024±97 


3967 


rj 


1196L46 


409 


792±69 


1052 




1678L74 


1898 


1050L113 


3991 


7r + 7T - 7T~ 


3654±199 


25208 


2300±187 


15723 


K*-K*° 


2030±98 


4878 


1298L130 


5672 


VP' 


4142±281 


20784 


2195L225 


17353 


Sum 


31302 ± 472 


70302 


18645±426 


63520 



first Gaussian. We require the invariant masses to be 
within ± 2.5cr (±2<r for the r\p~ mode) of Mu s - We 
have a total of 31302±472 tag candidates. Then we 
add a 7 candidate that satisfies our shower shape re- 
quirement. Regardless of whether or not the 7 forms 
a D* with the tag, for real D*D S events, the missing 
mass squared, MM* 2 , recoiling against the 7 and the 
D~ tag should peak at . We calculate 
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where Eqm (Pcm) ^ s ^he center-of-mass energy (mo- 
mentum), Eo 3 (jf D ) is the energy (momentum) of 
dividual widths and fi is the fractional area of the the fully reconstructed Dj tag, E 1 is the en- 
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ergy (momentum) of the additional 7. We use a 
kinematic fit that constrains the decay products of 
the D~ to Md s and conserves overall momentum 
and energy. All 7's in the event are used, except for 
those that are decay products of the Dj tag. 

The MM* 2 distribution from K+K~n~ tags is 
shown in Fig. O We fit all the modes individually 
to determine the number of tag events. This proce- 
dure is enhanced by having information on the shape 
of the signal function. We use fully reconstructed 
D~ D* + events, and examine the signal shape when 
one D s is ignored. The signal is fit to a Crystal Ball 
function [14| . which determines a and the shape of 
the tail. Though a varies somewhat between modes, 
the tail parameters don't change, since they depend 
on beam radiation and 7 energy resolution. 
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FIG. 2: The MM* 2 distribution from events with a 7 
in addition to the K + K~ix~~ tag. The curve is a fit to 
the Crystal Ball function and a 5th order Chebychev 
background function. 

Fits of MM* 2 in each mode when summed show 
18645±426 events within a ±2.5cr interval (see Ta- 
ble UJ. There is a small enhancement of (4.8 ±1.0)% 
in our ability to find tags in [i + i> (or tt + T>i>) events 
(tag bias) as compared with generic events. Addi- 
tional systematic errors are evaluated by changing 
the fitting range, using 4th and 6th order Chebychev 
background polynomials, and allowing the parame- 
ters of the tail of the fitting function to float, leading 
to an overall systematic uncertainty of 5%. 

Candidate fi + v events are required to have only a 
single additional track oppositely charged to the tag 
with an angle >35.9° with respect to the beam line. 
We also require that there not be any neutral en- 
ergy cluster detected of more than 300 MeV, which 
is especially useful to reject Df — > 7r + 7r° and 7771" + 
decays. Since here we are searching for events in 



which there is a single missing v, the missing mass 
squared, MM 2 , should peak at zero: 

MM 2 = (E CM - E Ds -E y - E^f (2) 

- (pcm - ~Pd s - % - %) , 

where E^ (p^) are the energy (momentum) of the 
candidate /i + track. 

We also make use of a set of kinematical con- 
straints and fit each event to two hypotheses: 

(1) the D~ tag is the daughter of a D*~ and 

(2) the D* + decays into jD+. The kinemati- 
cal constraints, in the center-of-mass frame, are 
~P~d s + ~&Di = 0' e cm = E Dg + Ed*, E d * = 
E C m/2+ (m 2 , - M 2 a ) /2E CM or E D . = E CM /2- 

(m 2 d , - M 2 s ) /2£ C M, M B% - M Ds = 143.6 MeV. 
In addition, we constrain the invariant mass of the 
Dj" tag to Md b - This gives a total of 7 constraints. 
The missing v four-vector needs to be determined, 
so we are left with a three-constraint fit. We perform 
an iterative fit minimizing \ 2 ■ To eliminate system- 
atic uncertainties that depend on understanding the 
absolute scale of the errors, we do not make a \ 2 cut 
but simply choose the 7 and the decay sequence in 
each event with the minimum \ 2 - 

We consider three separate cases: (i) the track de- 
posits < 300 MeV in the calorimeter, characteristic 
of a non-interacting pion or a fi + ; (ii) the track de- 
posits > 300 MeV in the calorimeter, characteristic 
of an interacting pion; or (iii) the track satisfies our 
electron selection criteria. The separation between 
muons and pions is not complete. Case (i) contains 
99% of the muons but also 60% of the pions, while 
case (ii) includes 1% of the muons and 40% of the 
pions 0] ■ Case (iii) does not include any signal but 
is used for background estimation. For cases (i) and 
(ii) we insist that the track not be identified as an 
electron or a kaon. Electron candidates have a match 
between the momentum measured in the tracking 
system and the energy deposited in the Csl calorime- 
ter, and dE/dx and RICH measurements consistent 
with this hypothesis. 

For the n + v final state the MM 2 distribution is 
modeled as the sum of two Gaussians centered at 
zero. A Monte Carlo (MC) simulation of the MM 2 
shows cr=0.025 GeV 2 after the fit. We check the 
resolution using the Df — > K K + mode. We search 
for events with at least one additional track identi- 
fied as a kaon using the RICH detector, in addition 
to a D~ tag. The MM 2 resolution is 0.025 GeV 2 in 
agreement with the simulation. 

In the 7r + 77j/ final state, the extra missing v re- 
sults in a smeared MM 2 distribution that is almost 
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triangular in shape starting near -0.05 GeV 2 , peak- 
ing near 0.10 GeV 2 , and ending at 0.75 GeV 2 . 
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FIG. 3: The MM 2 distributions from data using D~ tags, 
and one additional opposite-sign charged track and no 
extra energetic showers, for cases (i), (ii), and (iii). 

The MM 2 distributions from data are shown in 
Fig. [3] The overall signal region is -0.05 < MM 2 < 
0.20 GeV 2 . The upper limit is chosen to prevent 
background from r)ir + and K°ir + final states. The 
peak in Fig. [3]4) is due to D+ — > n + v. Below 0.20 
GeV 2 in both (i) and (ii) we have tt+vv events. 
The specific signal regions are: for n + v, —0.05 < 
MM 2 < 0.05 GeV 2 , corresponding to ±2<r; for tt+Vv, 
in case (i) 0.05 < MM 2 < 0.20 GeV 2 and in case (ii) 
-0.05 < MM 2 < 0.20 GeV 2 . In these regions we 
find 92, 31, and 25 events, respectively. 

We consider backgrounds from two sources: one 
from real D+ decays and the other from the back- 
ground under the single-tag signal peaks. For the 
latter, we estimate the background from data using 
side-bands of the invariant mass, shown in Fig. [TJ 
For case (i) we find 3.5 (properly normalized) back- 
ground events in the fi + v region and 2.5 back- 
grounds in the t+v region; for case (ii) we find 3 
events. Our total background estimate summing 
over all of these cases is 9.0±2.3 events. 

The background from real D+ decays is evaluated 
by identifying specific sources. For [i + v the only 
possible background is D + — > tt+tt . Using a 195 
pb _1 subsample of our data, we limit the branching 
fraction as < 1.1 x 10~ 3 at 90% C.L. 8]. This low 
rate coupled with the extra 7 veto yields a negligible 
contribution. The real Dj backgrounds for ir+Vv 



are listed in Table UTI Using the SM expected ratio of 
decay rates we calculate a contribution of 7.4 ix+vv 
events. 



TABLE II: Event backgrounds in the tt + w sample from 
real D + decays. 



Source 




B{%) 


case (i) 


case (ii) 


Sum 




Xfi+v 


8.2 


Q + l.S 

U -() 









jr+irV 


1.0 


0.03±0.04 


0.08±0.03 


0.11±0.04 




T+V 


6.4 








r+ 




1.5 


0.55±0.22 


0.64±0.24 


1.20±0.33 


r+ 


— » jl + UV 


1.0 


0.37±0.15 





0.37±0.15 


Sum 






i.oiF 


0.7±0.2 


1 7+ 1 - 8 

1 -0.4 



The event yield in the signal region, A^jot (92), is 
related to the number of tags, iV tag , the branching 
fractions, and the background A^bkgrd (3-5) as 



^Vdot - AWd = ^ta g ■ e[e'B(D+ -» n+v) 
+e"B(D+ -» n+Vv)], 



(3) 



where e (80.1%) includes the efficiencies (77.8%) for 
reconstructing the single charged track including fi- 
nal state radiation, (98.3)% for not having another 
unmatched cluster in the event with energy greater 
than 300 MeV, and the correction for the tag bias 
(4.8%); e' (91.4%) is the product of the 99.0% fj+ 
calorimeter efficiency and the 92.3% acceptance of 
the MM 2 cut of |MM 2 | < 0.05 GeV 2 ; e" (7.6%) is 
the fraction of tt + vv events contained in the /i + v sig- 
nal window (13.2%) times the 60% acceptance for a 
pion to deposit less than 300 MeV in the calorime- 
ter. Using B(t+ -> 7T+77) of (10.90±0.07)% 0, the 
ratio of the 7r + 77^ to fi + v widths is 1.059; we find: 

B(Df -> fi+v) = (0.594 ± 0.066 ± 0.031)%. (4) 



We can also sum the fj, + v and t+v contributions 
for -0.05 < MM 2 < 0.02 GeV 2 . Equation [3] still ap- 
plies. The number of signal and background events 
changes to 148 and 10.7, respectively, e' becomes 
96.2%, and e" increases to 45.2%. The effective 
branching fraction, assuming lepton universality, is 



2? eff (£>j 



(0.638 ±0.059 ±0.033)%. (5) 



The systematic errors on these branching fractions 
are dominated by the error on the number of tags 
(5%). Other errors include: (a) track finding (0.7%), 
determined from a detailed comparison of the sim- 
ulation with double tag events where one track is 
ignored; (b) the error due to the requirement that 
the charged track deposit no more than 300 MeV 
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in the calorimeter (1%), determined using two-body 
D° — » K~n+ decays [5j]; (c) the 7 veto efficiency 
(1%), determined by extrapolating measurements on 
fully reconstructed events. Systematic errors arising 
from the background estimates are negligible. The 
total systematic error for Eq. 0]is 5.2%, and is 5.1% 
for Eq. [3] as (b) doesn't apply here. 

We also analyze the t+v final state independently. 
For case (i) we define the signal region to be the in- 
terval 0.05<MM 2 <0.20 GeV 2 , while for case (ii) 
-0.05<MM 2 <0.20 GeV 2 . The upper limit on MM 2 
is chosen to avoid background from the tail of the 
K°n+ peak. The fractions of the MM 2 range ac- 
cepted are 32% and 45% for case (i) and (ii), respec- 
tively. 

We find 31 [25] events in the signal region with a 
background of 3.5 [5.1] events for case (i) [(ii)]. The 
branching fraction, averaging the two cases is 



We combine with our previous result f D + = 
222.6 ± 16. 7tji MeV @], and find 



B(D+ -> t+v) = (8.0 ± 1.3 ± 0.4)%, 



(6) 



where the systematic error includes a contribution 



of 0.06% from the uncertainty on £>(r H 



-77). 



We measure 13.4 ± 2.6 ± 0.2 for the ratio of t+v 
to /.i+v rates using Eq. [4] Here the systematic er- 
ror is dominated by the uncertainty on the mini- 
mum ionization cut. We also set an upper limit of 
B(D+ -> e+v) < 1.3 x 10~ 4 at 90% C.L. Both of 
these results are consistent with SM predictions and 
lepton universality. 

We perform an overall check of our procedures by 
measuring B(D+ -> K°K+). We compute the MM 2 
(Eq.[2]) using events with an additional charged track 
identified as a kaon. These track candidates have 
momenta of approximately 1 GeV/c; here the RICH 
has a pion to kaon fake rate of 1.1% with a kaon 
detection efficiency of 88.5% [12]. For this study, 
we do not veto events with extra charged tracks, or 
7's, because of the presence of the K°. We deter- 
mine B(D+ -> K°K+) = (2.90 ± 0.19 ± 0.18)%. 
This method gives a result in good agreement with 
preliminary CLEO-c results using double tags of 
(3.00 ± 0.19 ± 0.10)% 0; these results are not in- 
dependent. 

We also performed the entire analysis on a MC 
sample that is 4 times larger than the data sam- 
ple. The input branching fraction is 0.5% for fi+v 
and 6.57% for t+v, while our analysis measured 
(0.514±0.027)% for the case (i) [i+v signal and 
(0.521±0.024)% for y+v and t+v combined. 

Using B(D+ — > [i+v) from Eq. [5l and Eq. Q] with 
a D s lifetime of (500±7) x 10" 15 £ 



we extract 



f Dt /f D + = 1.23 ± 0.11 ± 0.04. 



(8) 



Lattice QCD predictions for f D + and the ratio 
f D + / fjj+ have been summarized by Onogi [16]. Our 
measurements are consistent with most calculations; 
examples are unquenched Lattice that predicts 249 ± 
3 ±16 MeV and 1.24 ± 0.01 ± 0.07 for the ratio 0, 
while a recent quenched prediction gives 266 ± 10 ± 
18 MeV and 1.13 ± 0.03 ± 0.05 [lj|. There is no 
evidence yet for any suppression in the ratio due to 
the presence of a virtual charged Higgs [Tol ] . 

The CLEO-c determination of f D + is the most 
accurate to date and consistent with other measure- 
ments @, H| . It also does not rely on the indepen- 
dent determination of any normalization mode (e.g. 
<fiir+). (We note that a preliminary CLEO-c result 
using D+ — » t+v, t+ — » e+vv [l9( is consistent with 
these results.) 
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f D + = 274 ± 13 ± 7 MeV. 



(7) 
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